Search for: All records

Abstract Applying full-waveform methods to image small-scale structures of geophysical interest buried within the Earth requires the computation of the seismic wavefield over large distances compared to the target wavelengths. This represents a considerable computational cost when using state-of-the-art numerical integration of the equations of motion in three-dimensional earth models. “Box Tomography” is a hybrid method that breaks up the wavefield computation into three parts, only one of which needs to be iterated for each model update, significantly saving computational time. To deploy this method in remote regions containing a fluid-solid boundary, one needs to construct artificial sources that confine the seismic wavefield within a small region that straddles this boundary. The difficulty arises from the need to combine the solid-fluid coupling with a hybrid numerical simulation in this region. Here, we report a reconciliation of different displacement potential expressions used for solving the acoustic wave equation and propose a unified framework for hybrid simulations. This represents a significant step towards applying ’Box Tomography’ in arbitrary regions inside the Earth, achieving a thousand-fold computational cost reduction compared to standard approaches without compromising accuracy. We also present examples of benchmarks of the hybrid simulations in the case of target regions at the ocean floor and the core-mantle boundary. 

SUMMARY We present a time-domain distributional finite-difference scheme based on the Lebedev staggered grid for the numerical simulation of wave propagation in acoustic and elastic media. The central aspect of the proposed method is the representation of the stresses and displacements with different sets of B-splines functions organized according to the staggered grid. The distributional finite-difference approach allows domain-decomposition, heterogeneity of the medium, curvilinear mesh, anisotropy, non-conformal interfaces, discontinuous grid and fluid–solid interfaces. Numerical examples show that the proposed scheme is suitable to model wave propagation through the Earth, where sharp interfaces separate large, relatively homogeneous layers. A few domains or elements are sufficient to represent the Earth’s internal structure without relying on advanced meshing techniques. We compare seismograms obtained with the proposed scheme and the spectral element method, and we show that our approach offers superior accuracy, reduced memory usage, and comparable efficiency. 

Abstract The existence of a thin, weak asthenospheric layer beneath Earth’s lithospheric plates is consistent with existing geological and geophysical constraints, including Pleistocene glacio-isostatic adjustment, modeling of gravity anomalies, studies of seismic anisotropy, and post-seismic rebound. Mantle convection models suggest that a pronounced weak zone beneath the upper thermal boundary layer (lithosphere) may be essential to the plate tectonic style of convection found on Earth. The asthenosphere is likely related to partial melting and the presence of water in the sub-lithospheric mantle, further implying that the long-term evolution of the Earth may be controlled by thermal regulation and volatile recycling that maintain a geotherm that approaches the wet mantle solidus at asthenospheric depths. 

ABSTRACT Seismic tomography is the most abundant source of information about the internal structure of the Earth at scales ranging from a few meters to thousands of kilometers. It constrains the properties of active volcanoes, earthquake fault zones, deep reservoirs and storage sites, glaciers and ice sheets, or the entire globe. It contributes to outstanding societal problems related to natural hazards, resource exploration, underground storage, and many more. The recent advances in seismic tomography are being translated to nondestructive testing, medical ultrasound, and helioseismology. Nearly 50 yr after its first successful applications, this article offers a snapshot of modern seismic tomography. Focused on major challenges and particularly promising research directions, it is intended to guide both Earth science professionals and early-career scientists. The individual contributions by the coauthors provide diverse perspectives on topics that may at first seem disconnected but are closely tied together by a few coherent threads: multiparameter inversion for properties related to dynamic processes, data quality, and geographic coverage, uncertainty quantification that is useful for geologic interpretation, new formulations of tomographic inverse problems that address concrete geologic questions more directly, and the presentation and quantitative comparison of tomographic models. It remains to be seen which of these problems will be considered solved, solved to some extent, or practically unsolvable over the next decade. 

Whether the two large low–shear velocity provinces (LLSVPs) at the base of Earth’s mantle are wide compact structures extending thousands of kilometers upward or bundles of distinct mantle plumes is the subject of debate. Full waveform shear wave tomography of the deep mantle beneath the Indian Ocean highlights the presence of several separate broad low-velocity conduits anchored at the core-mantle boundary in the eastern part of the African LLSVP, most clearly beneath La Réunion and Comores hot spots. The deep plumbing system beneath these hot spots may also include alternating vertical conduits and horizontal ponding zones, from 1000-km depth to the top of the asthenosphere, reminiscent of dyke and sills in crustal volcanic systems, albeit at a whole-mantle scale. 

Abstract Distributed acoustic sensing (DAS) is being explored in a variety of environments as a promising technology for the recording of seismic signals in dense array configurations. There is a particular interest for deploying DAS arrays on the ocean floor, presenting formidable challenges for conventional seismology. Taking advantage of the availability of a dark fiber on the Monterey Bay Accelerated Research System (MARS) 52 km offshore cable at Monterey Bay, California, in July 2022, we installed a DAS interrogator at the shore end of the cable with the intention of acquiring continuous data for a period of one year. Here, we describe the experiment and present examples of observations over the first six months of the deployment. 

Abstract We have extended the distributional finite difference method (DFDM) to simulate the seismic‐wave propagation in 3D regional earth models. DFDM shares similarities to the discontinuous finite element method on a global scale and to the finite difference method locally. Instead of using linear staggered finite‐difference operators, we employ DFDM operators based on B‐splines and a definition of derivatives in the sense of distributions, to obtain accurate spatial derivatives. The weighted residuals method used in DFDM's locally weak formalism of spatial derivatives accurately and naturally accounts for the free surface, curvilinear meshing, and solid‐fluid coupling, for which it only requires setting the shear modulus and the continuity condition to zero. The computational efficiency of DFDM is comparable to the spectral element method (SEM) due to the more accurate mass matrix and a new band‐diagonal mass factorization. Numerical examples show that the superior accuracy of the band‐diagonal mass and stiffness matrices in DFDM enables fewer points per wavelength, approaching the spectral limit, and compensates for the increased computational burden due to four Lebedev staggered grids. Specifically, DFDM needs 2.5 points per wavelength, compared to the five points per wavelength required in SEM for 0.5% waveform error in a homogeneous model. Notably, while maintaining the same accuracy in the solid domain, DFDM demonstrates superior accuracy in the fluid domain compared to SEM. To validate its accuracy and flexibility, we present various 3D benchmarks involving homogeneous and heterogeneous regional elastic models and solid‐fluid coupling in both Cartesian and spherical settings. 

Volcanic hotspots are thought to be fed by hot, active upwellings from the deep mantle, with excess temperatures ( T ex ) ~100° to 300°C higher than those of mid-ocean ridges. However, T ex estimates are limited in geographical coverage and often inconsistent for individual hotspots. We infer the temperature of oceanic hotspots and ridges simultaneously by converting seismic velocity to temperature. We show that while ~45% of plume-fed hotspots are hot ( T ex ≥ 155°C), ~15% are cold ( T ex ≤ 36°C) and ~40% are not hot enough to actively upwell (50°C ≤ T ex ≤ 136°C). Hot hotspots have an extremely high helium-3/helium-4 ratio and buoyancy flux, but cold hotspots do not. The latter may originate at upper mantle depths. Alternatively, the deep plumes that feed them may be entrained and cooled by small-scale convection. 

SUMMARY The presence of seismic anisotropy at the base of the Earth's mantle is well established, but there is no consensus on the deformation mechanisms in lower mantle minerals that could explain it. Strong anisotropy in magnesium post-perovskite (pPv) has been invoked, but different studies disagree on the dominant slip systems at play. Here, we aim to further constrain this by implementing the most recent results from atomistic models and high-pressure deformation experiments, coupled with a realistic composition and a 3-D geodynamic model, to compare the resulting deformation-induced anisotropy with seismic observations of the lowermost mantle. We account for forward and reverse phase transitions from bridgmanite (Pv) to pPv. We find that pPv with either dominant (001) or (010) slip can both explain the seismically observed anisotropy in colder regions where downwellings turn to horizontal flow, but only a model with dominant (001) slip matches seismic observations at the root of hotter large-scale upwellings. Allowing for partial melt does not change these conclusions, while it significantly increases the strength of anisotropy and reduces shear and compressional velocities at the base of upwellings.